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ABSTRACT 



Aims. Intrinsic alignments constitute the major astrophysical systematic for cosmological weak lensing surveys. We present a purely 
geometrical method with which one can study gravitational shear-intrinsic ellipticity correlations directly in weak lensing data. 
Methods. Linear combinations of second-order cosmic shear measures are constructed such that the intrinsic alignment signal is 
boosted while suppressing the contribution by gravitational lensing. We then assess the performance of a specific parametrisation 
of the weights entering these linear combinations for three representative survey models. Moreover a relation between this boosting 
technique and the intrinsic alignment removal via nulling is derived. 

Results. For future all-sky weak lensing surveys with photometric redshift information the boosting technique yields statistical errors 
on model parameters of intrinsic alignments whose order of magnitude is compatible with current constraints determined from indirect 
measurements. Parameter biases due to a residual cosmic shear signal are negligible in case of quasi-spectroscopic redshifts and 
remain sub-dominant for typical values of the photometric redshift scatter. We find good agreement between the performance of 
the intrinsic alignment removal based on the boosting technique and standard nulling methods, both reducing the cumulative signal- 
to-noise by about a factor of 6, which possibly indicates a fundamental limit in the separation of lensing and intrinsic alignment 
signals. 
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1. Introduction 



gravitational lensing: weak - large-scale structure of Universe - cosmological parameters - meth- 



Weak gravitational lensing of the large-scale structure is go- 
ing to be one of the major cosmological probes contributing 
to reveal th e properties of dark matter and dark ene rgy in the 
near future (lAlbrecht et alJ l2006: Pea cock et all I2006I) . Within 
the p ast decade the method has evolved from its first detec- 
tions ( iBacon et al.ll2000tlKaiser et al.ll2000t [van Waerbeke eTaTI 
2000; IWittman et al.ll2000l) to maturity, nowadays yielding sta- 
tistical constraints which a re compatible to other probes (for re- 
cent measurements see e .g. [Benjamin e t alj2007tlFu et alT2 008: 
Schrabback et al for a recent review see Munshi et al. 

2008). Planned surveys measuring weak lensing on cosmolog- 
ical scales, or cosmic shear in short, include Pan-STARRSQ, 
KIDfl DES0, LSSlQ and Euclirfl 

The increasingly large statistical power of these surveys de- 
mands a more and more thorough treatment of systematic errors. 
The major astrophysical contamination to cosmic shear is con- 
stituted by the intrinsic alignment of galaxies. To infer cosmic 
shear information from the correlation of galaxy ellipticities, it 
is usually assumed that the intrinsic shapes of galaxy images are 
purely random, so that only the desired correlations of gravita- 



tional shear (GG in the following) remain. However, due to in- 
teractions with the surrounding matter structure, galaxy shapes 
can intrinsically align, causing correlations between the intrin- 
sic ellipticities of galaxies (II hereafter). Moreover matter can 
influence the shape of a close-by galaxy via tidal forces and at 
the same time conttibute to the lensing signal of a background 
galaxy, thereby producing gravitational shear-intrinsic ellipticity 
correlations (GI hereafter). 

Intrinsic alignments have been subject to extensive stud- 
ies, b ot h analytical a n d usin g simulation s (ICroft & Metzlet 
2000; iHeavens et al.1 120001: iLee & Pen! 120001: iPenetal 
120001; ICatelan et al l 120011; ICrittenden et al.l bOOU Uind 120021; 
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van den Bosch et al. 2002t" Hirata & Seliakl 



Mack ev et al 

2004 iHevmans et all 120061: iBridle & Abdallal 12007 
Sembolon i et al.l 120081: lOkumura et al.l l2009t lOkumura & Jind 
2009; iBrainerd et al.ll2009l) . Results vary widely, but are mostly 
consistent with a contamination of the order 10 % by both II and 
GI signals for future weak lensing surveys, which can lead to 
serious biases on co smological parameters if left untreated (e.g. 
Bridle & King 2007). Intrinsic alignments depend intricately on 
the formation and evolution of galaxies within their dark matter 
environment, so that models cannot be expected to develop 
far beyond the current crude level in the near future. For the 
most recent advancement in intrinsic alignment modelling see 
ISchneider & Bridld(l2010l) . 

Using uncertain models of limited ac curacy for assessing 
systematics in statistical analyses is risky (Kitching et al. 2009). 
Therefore observational data which can put limits on the pos- 
sible range of intrinsic alignment signals are highly warranted. 
It should be noted that in principle intrinsic alignments consti- 
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tute an interesting cosmological signal worth investigating, shed- 
ding light onto the interaction between galaxies, their haloes, and 
the large-scale structure. Both II and GI cor relations have been 
subject to inves t igatio n s in several data sets (Brown e t al.ll2002t 



. 120041: [M andelba um et all ._ 
20071: iBrainerd et all l2009t iMandelb aum et alj 1201 Oh , results 



Hevmans et al 



. 200 6t iHirataetalJ 



ranging from null to significant detections, depending strongly 
on the type and colour of galaxies considered. 

However, none of these observations were direct measure- 
ments of intrinsic alignments for the galaxy populations and red- 
shifts which are most interesting for cosmic shear because in 
those cases the shear signal clearly dominates the correlations of 
galaxy ellipticities. While the II signal is observed at small red- 
shifts where cosmic shear is negligible, the GI term is usually 
inferred from cross-correlations between galaxy number densi- 
ties and ellipticities in samples with spectroscopic redshifts. The 
latter approach requires the assumption of a simple form of the 
galaxy bias, which is of limited accuracy and inapplicable on 
small scales. If one wishes to analyse larger galaxy samples for 
which only photometric redshift information is available, further 
signals such as galaxy- galaxy lensing contribute and need to be 
modelled carefully (see Bernstein 2009; Joachimi & Bridle 2010 
for an overview on the types of signals contributing to correla- 
tions between galaxy number density and ellipticity). 

The II signal is less of a concern because, in order to intrin- 
sically align, a pair of galaxies has to have interacted physically, 
and hence to be both close on the sky and in redshift. This fact 
can be used to remove II correlations dKing & Schneider! 2002, 
2003 ; iHevmans & Heavensl2003l iTakada & Whitd2004l) . partly 
in a fully model-independent way with only marginal loss of sta- 
tistical power if precise redshift information is available. The GI 
signal is not restricted to physically close pairs of galaxies, but 
it can also be eliminated in a purely geometr ical way via nulling 
techniques (Joachimi & Schneider 2008, 2009). However, a con- 
siderable loss of cosmological information is inherent to nulling, 
and hence, it is still desirable to have a reliable model of GI cor- 
relations at one's dis posal to be used with other methods control- 
ling this systematic dKin g 2005; B ridle & K ing 2007: lBernsteinl 
l2009tlZhang||20TollJoachimi & Bridldl2010h . 

In the following we will develop a model-independent tech- 
nique to extract the GI signal from a cosmic shear data set, 
thereby allowing for direct measurements of GI correlations on 
the most relevant galaxy samples. This 'GI boosting' approach 
can be regarded as complementary to nulling both in its purpose 
and in its implementation. Analogous to the nulling technique, 
we will construct linear combinations of second-order cosmic 
shear measures, making only use of the well-known characteris- 
tic redshift dependence of the GI and GG terms. 

This paper is organised as follows. In Sect.|2]we present the 
principle of GI boosting and derive general conditions, which 
are used in Sect. [3] to explicitly construct weight functions for 
the boosting transformation of the cosmic shear signal. Section 
[4] details the modelling which we apply in Sect. [5] to assess the 
performance of the boosting technique. In Sect.|6]we construct 
a method to remove GI correlations based on the GI boosting 
technique and investigate the relati on between the new approach 
and t he standard nulling method o f lJoachimi & Schneider! (T2008. 
2009), before we summarise and conclude in Sect.[7] 



2. Method 

2.1. Basic relations 

We will base our technique on a tomographic cosmic shear data 
set, i.e. correlations of galaxy ellipticities which are in addition 
split into subsamples according to the available redshift infor- 
mation. Analogous to the nulling technique the method outlined 
in the following does not affect angular scales, so that we can 
without loss of generality use tomographic power spectra as our 
two-point cosmic shear measures. For an ov erview on the basics 
of cosmic shear see e.g. ISchneiden (|2006) whose notation we 
mostly follow. 

The convergence power spectrum of cosmic shear, correlat- 
ing two galaxy samples i and j, reads 
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dX8 w (x)8 (J> (x) {l+z(x)YPs\-,X] . (1) 



where Pg is the three-dimensional matter power spectrum, I the 
angular frequency, and z the redshift. The integration runs over 
all comoving distances x U P t0 me comoving distance horizon 
X\m- Moreover we have introduced the lensing efficiency 
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where p^(x) is the probability distribution of comoving dis- 
tances for galaxy sample /. Note that we assume a spatially 
flat universe throughout. Similar to ([TJ, one can define a tomo- 
graphic pow er spectrum of shear-e llipticity correlations (for de- 
tails see e.g. Hirata & Seljak 2004]), 



2c 2 



J*™ dx {p ( Hx) g u \x) + g ( Hx) p (J \x)) 



{i+z{x))x~ l Psi\-,x 



(3) 



where P$i denotes the three-dimensional cross-power spectrum 
between matter density contrast and intrinsic shear fielcQ. Only 
one of the terms in (0 is non-vanishing unless the probability 
distributions overlap. As II correlations can readily be removed 
before applying a treatment of the GI signal, we neglect them 
in this work, so that the total power spectrum, i.e. the actual ob- 
servable in our study, is given by 



(4) 



A discussion on how II correlations affect the boosting technique 
is provided in Sect. [7] 

To derive expressions for the transformed signals, we as- 
sume that precise redshift, or equivalently distance, information 
is available, so that the survey can be sliced into thin tomo- 
graphic bins. One can then approximate p^(x) ~ 8d(x ~ Xi)> 
where Xi is an appropriately chosen comoving distance in bin i. 
Here <5d denotes the Dirac delta distribution. The lensing effi- 
ciency (f2]l can then be written in the form 



(Xi) -» g(Xj>Xi) = 



1-f if Xi < Xi 







else . 



(5) 



6 The intrinsic shear is defined as the correlated part of the 
intrinsic ellipticity of a galaxy image (e.g. iHirata & Seliakl 120041 : 
Joachimi & Bridle 2010). One can then proceed to construct an intrin- 
sic shear field by assigning to every point in space the intrinsic shear a 
galaxy would h ave at this positio n. For instance, if the intrinsic align- 
ment model of Catelan et al. (2001) held true, this could simply be done 
by computing the quadrupole of the local gravitational field. 
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With these approximations the power spectra (0 and (0 turn 
into 



9H*Ql f™ n te*> 



4c 4 



dx g(Xi,x) g(Xj,x) (6) 



x {i+z(x)rp s ^-,x\ ; 



3flgQ m ( l+zfo) n 

^GiCr/,^'^) = —^—{giXjiXi) — — — Pa\—,Xi\ ( ? ) 

+ giXhXj) Psi —,Xj } > 



where the dependence of the power spectra on the comoving 
distances of the two galaxy samples involved was made ex- 
plicit. Note that if Xt < Xh on ty tne mst term contributes to 
Pgi(X uXji whereas for^, > Xj on ly the second term is non- 
zero. 



2.2. Signal transformation 

We seek to find linear combinations of tomographic second- 
order cosmic shear measures such that in the resulting measures 
the cosmic shear signal is largely suppressed with respect to the 
GI signal. The starti ng point is analogous to the n ulling tech- 
nique as outlined by Joachimi & Sch neider! d2008l) . We define 
transformed power spectra as 



55 L ' 



d X B & (x) Pobs(Xi,X,-0 



(8) 



where B ( '^(x) is a weight function yet to be determined. Note that 
([8]) holds also for both the GG and GI contributions individually 
as the observed power spectrum is a linear superposition of the 
two, see (0. We will investigate two different choices for the 
lower boundary of the integration Xmm in this work. To construct 
the boosting technique, we choose the maximum range Xmm - 
whereas in Sect. l6.ll we will set^ m ; n = X i instead. 
Inserting (0 into the definition (0, one finds that 



<(0 = 



3H 2 Q„ 

2c 2 



dx B (i) (x) \g(x,Xi) 
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where we defined the function 



G (i) (x) 



= dxB'Hx)^-^ 



(10) 



Note that the integration absorbed into G^'Hx) starts at,y, which 
corresponds to a lower boundary of X i in the integral over the 
first term in (0. This can be done because g(x,xd vanishes for 
X < Xi' see 0. Likewise, x < Xi holds for the second term in 
(0 to be non-zero, so that the upper boundary of this integration 
is changed to Xi- In addition 0, with only its first term non- 
vanishing, can be inserted. The first term in the final expression 



of (0 is generated by GI correlations originating from matter 
at the distance^,. Note that in the approximation of thin redshift 
slices which we are working in this term is just a rescaled version 
of 0. Due to our choice Xmm - 0, the transformed GI signal 
receives a further contribution from shear-ellipticity correlations 
generated ?&x < Xu collected into the second term of (0. 

Transforming the lensing signal analogously by plugging (0 
into 01, one arrives at 



n&(0 = 



4c 4 
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X, 



4*-|i- — \G (i) (x){i+z(x)} 2 P s \-,x) 



\X 



Again, < TT~Qb was used to produce the final expression. The con- 
ditions x < Xi an d X > X< imposed by 0, result in the upper 
boundary of the first and the lower boundary of the second in- 
tegral, respectively. The transformed cosmic shear signal thus 
depends on the form of G (l) (x) in the interval between and^,. 
To suppress the GG signal, G^(x) should be chosen such that 
the integral in the final expression of (fTTT i is close to zero while 
at the same time G (,) 0t,) has to be comparatively large to boost 
the GI contribution, see (0. 

In reality line-of-sight information will not be available in 
terms of comoving distances, but rather in terms of the observ- 
able redshift. Furthermore the galaxy redshift distributions will 
have a finite width and also overlap due to scatter, in particular 
if only photometric redshift information is available as will be 
the case for the vast majority of galaxies in future cosmic shear 
surveys. To arrive at a practical prescription for constructing the 
transformed power spectra, we therefore change the integration 
variable in (0 to redshift and subsequently discretise the inte- 
gral, yielding 



£ B (i \x(zj))P^)x'(zj)Azj 



(12) 



where x'(z) is the derivative of comoving distance with respect to 
redshift, and Azj is the width of redshift bin j. In total galaxy 
samples are available for study. Here and in the following we 
identify Zi = z(xd- The condition Xmm = used for the boosting 
technique translates into j m [ n = 1 . 

2.3. Solving for the weight function 

In the foregoing section we saw that the GI signal can be 
boosted, and the GG signal at the same time suppressed, by 
formulating conditions on the function G^(x). Via its defining 
equation ( fTUb it is related to the weight function B^(x) that en- 
ters the transformation (0. Hence, to obtain a boosting trans- 
formation, one has to solve ( TT0T > for B w 0f) for a given function 
G (i) (x). 

We begin by noting that ( fTOb is a Volterra integral equation 
of the first kind. It has a kernel that is linear in the integration 
variable, so that one can readily solve for the weight function by 
differentiating twice, resulting in 

( , d 2 G®(Y) 
B '(x) = X 



dx 2 



(13) 



We have found the solution of the inhomogeneous Volterra equa- 
tion ( [Tot under the premises that G ij \x) is twice continuously 
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differentiable, G w 0fhor) = and dG®/dx\ Xbal = 0. If one spec- 
ifies G w 0f) down to a value Xmm, then B®(x) is well-defined 
in the range [^ m i„,^hor] by ( fTOb . Note that if we dropped the 
assumption of a flat universe, ( fTOb would still be solvable, but 
analytical progress would be hampered. 

To find the solution of the homogeneous equation, obtained 
from ( fTOb by setting G (,) (y) = 0, we define 



b(x) = — ^ H( Xhm -x) ; fix) = XH(X) , 



X 



(14) 



where H(x) denotes the Heaviside step function. Then ( fTOb can 
be re-written as a cross-correlation, 



G & (X) 



f 



dxb(x)f(X-X) = {b*f) (X) 



(15) 



The introduction of the Heaviside functions in (fT4l was used to 
extend the integration to zero and infinity. If we denote Fourier 
transforms by a tilde, the convolution theorem yields G w = b f. 
From this equation is it readily seen that for G (l) (x) = it follows 
B^ix) = in the interval [^,^hor]- Hence the solution of the ho- 
mogeneous Volterra equation consists only of the trivial one and 
( fT3T > constitutes the full, unique solution of ( fTOb . In summary, for 
a given G (,) (x) that fulfils the conditions imposed by (|9]l and (fTTT) . 
we can calculate the corresponding weight function via (TT~3T > and 
use the result to construct transformed power spectra (8). 

Note the analogy between ( fTOb and the definition of the lens- 
ing efficiency (|2]). This can be interpreted as G^'fx) being a mod- 
ified lensing efficiency, which is then used to construct an alter- 
native lensing convergence with desired properties chosen via 
G (l) (x)- For details on this view see the motivatio n of the nulling 
technique given in I Joachimi & Schneider! d2008h . 



3. Construction of weights 

Apart from the requirements formulated in Sect. 12. 21 to ensure a 
boosting of the GI signal with respect to cosmic shear, the choice 
of G w 0t) is arbitrary. In the following we choose a specific 
parametrisation of G (l> (x) which is convenient and intuitive, but 
not necessarily optimal. Its base is a Gaussian that is peaked at 
Xi, which fosters a strong contribution of GI correlations via the 
first term of (0. Some additional flexibility is needed atx < %u 
allowing for sign changes of G w Ct) to downweight the lensing 
signal. We define 



G®(x) = N ex?{- ^ f m)2 ^ v -/M , 



(16) 



where N, cr, b, and^f m are free parameters. All four parameters 
depend on the choice of galaxy sample i, but we do not spec- 
ify this dependence for reasons of better readability. The first 
derivative of G^(x) with respect to comoving distance reads 



dG (i ) 
dx 



■(x) = N exp 



(X-Xm) 2 



[l-2(x- 



b)^^ \ . (17) 



From this result and by means of ( fT3l l one readily obtains the 
weight function 



B (i \x) = N^r exp^- 



2x 



(X-XmY 



(18) 



x {2(x-b) (X \ S -3x + 2 Xm + b 




Fig. 1. Parameter b as a function of cr, i.e. the width of the 
Gaussian in G^(x). Plotted are the results for different peak po- 
sitions Xi as indicated in the legend. For most of the considered 
range, (l23b provides an excellent approximation. Only for small 
Xi in combination with large cr do the curves start to level off. 
Note that b, cr, and^-,- are given in abstract units of 1 in this plot. 

The normalisation of G^ix) is related to the one of B^(x) via 
([Tol l, but is otherwise irrelevant to the problem. We fix N by 
requiring 

N 

I dx {B (i) (x)} * J] {B (i) (x(Zj))} x'(Zj) Azj = 1 . (19) 

Note that since N depends on the other free parameters, e.g. cr, 
a consistent normalisation is actually important when studying 
G w 0f) as a function of these parameters, as we will do in Sects. [5] 
and0 

Two of the remaining three free parameters of ( [ToT i will now 
be used to boost (0 and suppress ( fTTT i. First, we demand that 
(fT6l l is peaked at^-,-, i.e. dG^/dx \ Xi = 0. Using (JTTJ, we obtain 



Xm =Xi 



(Xi - by 



(20) 



The second condition should render the integral in (fTTT) close to 
zero. While it is possible to numerically determine for instance 
the parameter b such that this condition is fulfilled for every an- 
gular frequency individually, we prefer to proceed in a way that 
does not rely on a model of cosmic shear power spectra at all. 
We note that if the width of the Gaussian cr is relatively small, 
the support of the integral in (fTTT) has a small range and hence 
Pg can be well approximated as only varying slowly. The depen- 
dence on redshift should be roue hly P s (k,z) oc D(z) 2 oc (1 + Z )- 2 , 
where D(z) is the linear growth factor for which we assumed 
D(z) oc (1 + z) _1 as holds true in the matter-dominated epoch. 
This redshift dependence then cancels the (1 + z) 2 term in (fTTT i. 
so that we consider the condition 



G W ( X ) = . 



(21) 



Inserting ( fToT i together with (1201 1. and making the further defi- 
nitions y = (x - Xi)/°~ an d m = ix i ~ b)/cr, we transform this 
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The approximation in the first equality refers to replacing the 
lower boundary of the integral ~Xilo~ by — oo, which is valid if 
Xi/cr » 1, i.e. if G^ix) is compact (cr «: 1) and peaks not too 
close to zero (xi » 0). The root of the term in curly brackets can 
be found numerically, resulting in 



b(cr,Xi) — Xi~ mo ~ w i tn m ~ 1-10687 



(23) 



0.9 



1.1 



1.2 



We have solved (f2TT > directly and plot the resulting b in Fig.Q] 
We find excellent agreement with the approximate solution (1231 
as long as the assumption discussed above is fulfilled. Significant 
deviations from the linear behaviour of b as a function of cr 
are only found for Xi/°~ ^ 2. For reasons of simplicity we will 
restrict ourselves to cases where the approximation (1231 holds. 
This means in particular that we will not consider signals at very 
small redshifts, where Xi is necessarily small. In practice we use 
the condition G w (^(z m m)) ~ with z m in the minimum redshift 
used in the survey as a simple cross-check to ensure that this 
approximation is sufficiently accurate. 

The conditions specified above are strictly fulfilled only for 
continuous x or Z- However, we will in practice use the dis- 
cretised transformation (fT2l and thus have to make sure that 
GI boosting and GG suppression work accurately also in this 
case. Via a procedure outlined in the following, we optimise 
the remaining free parameter cr to guarantee a good sampling of 
G^(x) by the discrete set of weights B^HxiZj)) with j = 1 , .. , N z , 
thereby fulfilling dG®/dx \ Xi - and (1211 to good accuracy. 

As the sampling points of ( fT2l we choose the medians of 
the redshift distributions of the galaxy samples employed. It is 
expected that the optimal choice of the parameter cr, denoted 
by cr op t in the following, will depend intricately on the posi- 
tions of these sampling points and hence on the redshift distri- 
butions of the different galaxy samples in the cosmic shear data, 
in particular if the number of sampling points is small, e.g. if 
the distributions have a large scatter. Since the binning is done 
in terms of redshift, it is convenient to work with the quantity 
cr- = {x'(Zi)} 1 cr instead of cr. We will also give our choices of 
cr opt in terms of cr z throughout. 

We introduce the discrete version of the function G®(x), 



Fig. 2. Functions G w (y(z)) and B^(x(z)) for different Zu as in- 
dicated by the vertical grey lines, and redshift binnings, using 
the width o~ z determined via (1251 in each case. In addition the 
sampling points corresponding to the median redshifts of the 
bins are shown for B il \xiz)). The normalisation has been cho- 
sen according to d!9t . See Table[T]for an overview on the survey 
models referred to in the following. Top panel: For the spectro- 
scopic survey S and Zi = 0.53; cr. = 0.055. Centre panel: For 
the survey with good photometric redshifts PI and z, = 0.76; 
cr z = 0.085. Bottom panel: For the survey with standard photo- 
metric redshifts P2 and z, = 0.98; cr z = 0.055. 



integral as follows, 



\G {1) (X) 



o \ Xi 
No 3 <*0 



(22) 



Xi 
No- 3 
8m 2 Xi 



f 

•J — o 



dyy(y + m) ex Pi~\y + ^\ 



2m (1 - 2m z ) e 



2x .-l/(4m) 2 



1+Erfl — 

Lm 



G' (i) (x(z k )) 



= ^ B U(x(Zj))(l~^x'(Zj)Azj. (24) 



Then we consider the root mean square deviation of all function 
values G'®(x(zkJ) used, 



1 Nz 

r{cr z ) = . — ^ \G'®(x(z k ), cr z ) - G(%(Z*), cr z 
* z k=l 



(25) 



as a criterion for how well G w (y) is sampled by the discrete set 
of function values B®(x(ZjJ) entering d24l i. In the equation above 
we have made the dependence on cr z explicit in the arguments. 
We emphasise that the determination of cr z via the diagnostic £ is 
optimal only in the sense that it allows us to find a representative 
sampling of G {i) (x) such that 8G®/dx\ Xl = and (gT) hold to 
good accuracy. It will in general not yield an optimal amplifica- 
tion of the GI signal over the lensing signal, which depends on 
the explicit form of both signals. 

In Fig. [2] we have plotted a selection of typical results for 
G w (x> and the corresponding weight function B (/) 0f). As com- 
mon features of G^(x) a distinct peak at z, and a negative dip 
at z < Zi, the latter necessary to fulfil (|2TT i. are discernible. The 
weight function B^(x) has three pronounced extrema of which 
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the central one is located at Zu plus a shallow fourth one at low 
redshift. 

Note that the method laid out here is completely independent 
of any assumptions about the angular dependence of both the 
underlying lensing and intrinsic alignment signals. To determine 
the weights entering ( fI2t . we only make use of the well-known 
redshift dependence of the GI and GG signals, plus the redshift 
binning of the survey to be analysed. We note that the weights 
B^(x(Zj)) depend on £2 m and possibly further cosmological pa- 
rameters via the distance-redshift relation. However, the same 
applies to the weights used i n the standard nulling techn ique, 
and from the investigation by Joachim i & Schneiderl d2009l) we 
conclude that this dependence is weak and that the assumption 
of an incorrect cosmology when constructing the weights is un- 
critical. 



4. Modelling 

To assess the performance of the boosting technique, we need 
to model both the cosmic shear and the intrinsic alignment sig- 
nals. To this end, we assume a spatially fiat ACDM universe 
with matter density parameter Q m = 0.25 and Hubble parame- 
ter h = 0.7. The matter power spectrum has a primordial slope 
n s = 1 .0 and normalisa tion erg = 0.8. The tra nsfer function is 
computed according to lEisenstein & Hvj (1 19981) . using a baryon 
density parameter of £2b = 0.05, while the non-linear evolu- 
tion of the power spectrum is determined by the fit formula of 



ISmith et 



le power i 
alj(|2003) 



We use the linea r alignment model (ICatelan et alJ 1200 it 
iHirata &~S eliak 200% to calculate the matter density-intrinsic 
power spectrum, 



Pgi (k, z) = -C G i p c 



n m (i+z) 2 

D(z) 



PsiKz) 



(26) 



where the normalisation is c hosen such that CaPcr ~ 0.0134 
dJoachimi & Schneiderl 12009, and references therein) with p cr 
the critical density. In (l26b we use the full non-linea r matter 
power spectrum as suggested by Bridle &Kingl (120071) . While 
this conjecture lacks a s ound physical basis, the resulting sig- 
nal fits existing data well dBridle & Kingll2007b and has recently 
been shown to also be co nsistent with halo model calculations 
(ISchneider & Bridldl20loh . 

A cosmic shear survey is modelled b y assuming an ov erall 
galaxy redshift distribution according to Sma il et aTT (fr994). 



(27) 



with zo = 0.64 and (3 — 1.5, corresponding to a median redshift 
of Zmed = 0.9. We cut the distribution below z m in =0.2 and above 
Zmax =2.0 and normalise ( l27l i in that interval. The overall red- 
shift distribution is then sliced into disjoint bins. In those cases 
where a scatter due to photometric redshift estimates is present, 
we assume the distribution of photometric redshifts for a given 
true redshift to be a Gaussian, centred on the true redshift and 
with a width of <x ph (l + z). The distributions of true redshifts 
/? w (z) for each photometric red shift bin i are then computed ac- 
cording to a scheme detailed in Joachi mi & Schneiderl (120091) . 

We consider three different survey models which are sum- 
marised in Table Q] All of these surveys are assumed to cover 
the whole extragalactic sky, i.e. A survey = 20, OOOdeg 2 . To cal- 
culate shape noise, we use an intrinsic ellipticity dispersion of 
cr e = 0.35 throughout. 



Table 1. Overview on the different survey models used. 
identifier redshifts bin width tr ph n g [arcmin -2 ] 



S 
PI 
P2 



spectroscopic 
good photo-z 
standard photo-z 



0.01(1 +z) 
0.01(1 +z) 
0.02(1+*) 




0.03 
0.05 



1 

10 
40 



First, we construct a 'spectroscopic' survey S for which red- 
shift bins are assigned with width 0.01(1 + z) and no scatter. 
In this case the signals are calculated to excellent approxima- 
tion not over the complete bin width, but at the median red- 
shifts of each bin. Whilst it is in principle possible to achieve 
such a dense redshift binning and small scatter with photomet- 
ric redshifts (see e.g. Illbert et all 120091) . it is more likely that 
future large-area spectroscopic surveys fit into this category. In 
any case the number of available galaxies will be small. Taking 
the wide spectroscop ic survey of the Euclid mission as reference 
dLaureiis et al.ll2Q09h . we set the overall galaxy number density 
to « g = 1 arcmin -2 . 

Second, we create a survey that features high-quality pho- 
tometric redshift data, termed PI. We choose the same binning 
scheme as for the first case, but introduce a photometric red- 
shift scatter of <x p h = 0.03, corresponding to the target value of 
the Euclid imaging survey. To be conservative, we assume that 
this photometric redshift quality is only attainable for a subset 
of galaxies and set n g — 10 arcmin -2 . Finally, we make use of a 
setup P2 with redshift binning in steps of 0.02(1 + z) and scatter 
crph = 0.05, which can be regarded as representative of a stan- 
dard future imaging surve y designed to do cosmic shear. Again 



referring to lLaureiis et alJ d2009h . we adopt n g = 40 arcmin 2 in 



this case. 

The photometric redshift bin widths are chosen such that the 
associated distributions of neighbouring bins can still be well 
distinguished. We have found that narrowing the bin widths sub- 
stantially below about l/3<x p h deteriorates the performance of 
the boosting technique. It should be noted that spectroscopic red- 
shifts as well as photometric redshifts of high quality are usually 
limited to a brighter subset of galaxies, therefore altering the 
overall redshift distribution of galaxies. However, to facilitate 
the comparison between the three survey models under scrutiny, 
we keep piot(z) as specified above. 

With the three-dimensional GG and GI power spectra and the 
redshift distributions p^'fx) = p <l \z)/x'(z) at hand, one can cal- 
culate the tomographic power spectra according to (Q]) and (0. 
For the further analysis we divide the angular frequency range 
into N t = 200 logarithmic bins between £ = 10 and I = 20000. 



5. Performance of GI boosting 

5.1. Boosted signals 

To condense the performance of the boosting technique into a 
single number, we define the median with respect to angular fre- 
quency of the ratio of GI over GG signal, 



rQi = median 



(28) 



where X can be replaced by any tomography power spectrum 
p( l fi(€) or the transformed power spectra n w (f). Note that this 
quantity is not available from a real survey because we are not 
able to separate the GG and GI signals, but only extract their sum 
from the data. We have chosen the median in ( l28l since we find 
that the mean is not a robust measure for two reasons. First, if the 



B. Joachimi and P. Schneider: Intrinsic alignment boosting 



7 




0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 
a. 

Z 

Fig. 3. Diagnostic see d25l ). and GI over GG ratio ra, see d28l ). 
as a function of cr z . Shown is rci as solid curve and f as dashed 
curve. Since the normalisation of £ is arbitrary, we have rescaled 
£ for easier inspection. In each panel the choice of cr opt is marked 
with an arrow. Note that this choice was made without resorting 
to rci which is not measurable from real data. Upper panel: For 
the spectroscopic survey S at z, = 0.53. In addition we have 
plotted rci for the case <x p h = 0.03. Note that ( remains the same 
in both cases. Centre panel: For the survey PI (cr p h = 0.03) at 
Zi = 0.76. Note the dip of £ at <r z ~ 0.04 which is caused by 
the weight function B^(x) being sampled close to its extrema. 
We have added the curves for zi = 0.81 as grey lines, where 
B^(x) is sampled almost exactly at its extrema, leading to a local 
maximum in rci at o~ z = 0.038, traced well by a corresponding 
minimum in (. Lower panel: For the survey P2 (cr ph = 0.05) at 
Zi = 0.98. 



GG signal is suppressed by several orders of magnitude, numer- 
ical noise stemming from the computation of the power spectra 
can become important, leading to unphysical dips in the residual 
power spectrum. Second, the residual GG signal may have sign 
changes close to which ra becomes very large, thus dominating 
the mean. Both effects would mimic a stronger boosting than is 
actually observed. 



Table 2. Summary of rci for different values of z, and the three 
survey models. Given are values of rci of the original power 
spectra for the auto-correlation ( 'auto' ), the cross-correlation be- 
tween bin i and the background bin with index (i+N z )/2 ('mid'), 
and the cross-correlation between bin i and the most distant bin 
with index N z at z < 2 ('far'). The tag 'boost' stands for the 
transformed signals. In addition <x opt is listed for every consid- 
ered case. 



survey 


Zi 


O" opt 


r G i(auto) 


r G i(mid) 


rci (far) 


r G i (boost) 


S 


0.53 


0.055 


0.04 


0.47 


0.56 


78.26 




0.76 


0.085 


0.02 


0.25 


0.32 


415.79 




0.96 


0.095 


0.01 


0.16 


0.22 


62.10 


PI 


0.53 


0.055 


0.12 


0.46 


0.54 


13.56 




0.76 


0.085 


0.06 


0.24 


0.31 


20.49 




0.96 


0.095 


0.04 


0.16 


0.21 


25.86 


P2 


0.52 


0.045 


0.19 


0.47 


0.55 


5.97 




0.74 


0.050 


0.10 


0.26 


0.32 


7.75 




0.98 


0.055 


0.06 


0.15 


0.20 


11.63 



In Fig. [3] we show rci, together with the diagnostic ( as de- 
fined in S251 . as a function of cr z for one z,- per survey model. 
Overall we find that small values of ( indeed indicate regimes 
of cr z in which the GI signal is well boosted. It is important to 
note that the absolute value of ( is meaningless due to the ar- 
bitrariness in the overall amplitude of G®(y). When G w 0f) is 
no longer well sampled for small <x-, £ features a clear increase. 
Sometimes secondary minima in £ can be observed, see the cen- 
tre panel of Fig. [3] which is caused by the sampling points being 
consecutively placed at the extrema of B^(x). Thereby, although 
only sparsely sampled, the discrete form ( l24"j i captures the main 
characteristics of B^'Xx) an d hence can well represent G w 0t0, 
yielding a small value of ( . 

In the top panel of Fig. [3] rci for both surveys S and PI is 
given. Since the binning scheme is identical for both surveys, ( 
is the same. This example demonstrates that rci depends consid- 
erably on the details of the actual signals, in this case a change 
from <r ph = to cr ph = 0.03. The diagnostic ( does not trace 
the boosting of the actual signals and can consequently not be 
exploited to find the maximum rci- However, for both surveys ( 
identifies the regime of small cr z in which the boosting performs 
worse and which thus should be avoided. In the case <x p h = 0.05 
the sampling in redshift becomes fully insufficient for small cr z . 
Accordingly, ( rises sharply, and the GG signal starts to domi- 
nate again. 

The optimal width of G w 0t) can be chosen freely in the in- 
terval where £ is stable and small. If there is a clear minimum, 
we place cr opt there; otherwise we set cr opt to a small value in the 
interval where f is small, see e.g. the centre panel of Fig. [3] This 
assignment of cr opt may not be unique, but it is uncritical. Note 
that the weight functions corresponding to the optimum cases of 
the examples shown in Fig. [3] are those depicted in Fig. [2] We 
emphasise again that rci cannot be measured from real data, and 
accordingly we do not use this quantity to determine <x opt . 

One might expect that the denser the sampling points of 
G w 0f) and B^(x) can be placed, the more sharply peaked weight 
functions can be well represented by the discrete sampling, and 
thus smaller values of <t z could be chosen. However, consider 
the case <x ph = 0.05 and z, = 0.98 which is shown in the bottom 
panels of both Fig s. [2] and [3] Although cr opt is small compared 
to e.g. our findings for survey PI, the sparse sampling obviously 
captures the main features of the weight function and hence re- 
sults in a small (. 
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Fig. 4. Left column: Example set of lensing (GG) and intrinsic alignment (GI) tomography power spectra for the spectroscopic 
survey S and zi = 0.53. The GG signal is shown as solid line, the GI signal as dotted line. The upper three panels show power 
spectra for different background bins /', i.e. auto-correlations (j = i, 'auto'), cross-correlations with a bin at intermediate redshift 
(j = (i + Ay/2, 'mid'), and cross-correlations with the most distant bin (j = N z , 'far'). In the bottom panel the transformed GG 
and GI signals are plotted ('boost'). Note that absolute values of the power spectra are shown throughout. Centre column: Same 
as above, but for the survey PI (cr ph = 0.03) and Zi = 0.76. Right column: Same as above, but for the survey P2 (cr p h = 0.05) and 
Zi = 0.98. 



In Fig. |4] we have plotted example sets of original and trans- 
formed power spectra, the latter each computed for the optimal 
values of <x~. Table [2] lists the corresponding values of ra and 
cr pt for the three survey models and three redshifts z; each, in- 
cluding the cases depicted in the figure. The GI over GG ratio 
rci for the original power spectra ranges from about 1 % to 50 %. 
For a correlation between galaxy samples i and j with z, < Zj, 
ra increases strongly with the separation between z / and zu Both 
the GI and GG signals show this behaviour due to y]). Since the 



cosmic shear signal is generated by all the matter between z — 
and Zi with the highest efficiency at z,-/2, whereas the intrinsic 
alignment contribution stems form matter around zu me GI sig- 
nal has the stronger dependence on redshift, causing the increase 
in ra- For the non-linear version of the linear alignment model 
this effect can lead to a GI signal whose absolute value can come 
close to or e ven surpass the cosmi c shear signal for large Zj - Zi, 
see e.g. also lBridle & Kingl(l2007h . 
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Table 3. Statistical errors <x stat and residual biases b sys for the dif- 
ferent survey models used. The left column shows marginalised 
1 cr errors on the amplitude A of the GI signal and on y, quan- 
tifying an additional redshift dependence. In the right column 
results obtained for only varying A are listed. 



survey 


parameter 


y varied 


y 


fixed 






0"stat 


o sys 


0"stal 


b sys 


S 


A 


2.885 


-0.009 


0.827 


-0.004 




1 


7.356 


0.012 






PI 


A 


0.712 


-0.081 


0.172 


-0.046 




7 


1.776 


0.090 






P2 


A 


0.697 


-0.181 


0.272 


-0.171 




1 


2.017 


0.031 







It is evident from Table [2] that the better resolved the red- 
shift information is, the more can the GI signal be boosted. For 
quasi-spectroscopic data, the residual GG contribution is well 
below the 2 %-level and hence expected to be negligible. In the 
case Zi = 0.76 we find by chance a near-total cancellation of the 
cosmic shear signal. For good photo-z data with <x ph = 0.03 the 
method is also effective, yielding rci well in excess of 10, so 
that any biases due to the residual GG contribution are likely to 
remain below the statistical errors of intrinsic alignment param- 
eters. For survey P2 it is still possible to produce a dominating 
GI signal, with ra between approximately 6 and 12, but a GG 
residual exceeding 10 % may require further treatment to avoid 
a significant bias. 

5.2. Parameter constraints 

The boosted GI signal has the potential use of directly con- 
straining models of intrinsic alignments, provided that the sta- 
tistical power is sufficiently high and that systematics due to 
residual GG contributions are under control. We set up a simple 
intrinsic alignment m odel and use the Fisher matrix formalism 
dTegmark et al.|[l997h to forecast expected errors and biases on 
its free parameters. We define 



^model 



{k,z)=AP sl {k,z) 



1 +z 

1 + Znix 



(29) 



where Psi is given by the linear alignment model (f26b . The free 
parameters are A and y, i.e. we allow for an arbitrary signal 
amplitude and an additional redshift dependence. The fiducial 
model is (l26l l, so A — 1 and y = 0, and we set Zmv _ 0.3. 
The s ame parametrisation was e.g. used by iMan delbaum et al.l 
(120101) . 

Assuming that the signal covariance is itself not parameter- 
dependent, which holds to very go od accuracy for the large sur- 
vey we consider dEifler et al.l l2009). the Fisher matrix reads 



p 



zz 

t u 



dp? 



cov- 1 (n®(0,iig(/)) 



(30) 



for a parameter vector p = {A,y}. Using (fT2l . one can readily 
relate the covariance of the transformed power spectra to that of 
the original power spectra, 



cov (n® (o, n$(0) = J] B(r) 0dzk)) (31) 

k,t=0 

x Cov P ( i{t)) X '(z k )x'(zi) Azu Az, . 



The power spectru m covariance in turn is given by (see 
[ et al.ll2008l ; 



Joachimi i 



Cov(p^), 



and references therein) 



GI 



{€)) = 



2n 



A t\t 

-^survey L L -^ 1 



GI 



(32) 



+ p^({)p ( £\o\ withp^ = p; 



GI 



,07) 



aft/) 
GI 



°1 



where A£ is the width of the angular frequency bins and n w the 
number of galaxies belonging to sample i. Equation d32b holds 
under the assumptions of Gaussian density fluctuations, a uni- 
form sampling of galaxies, and a simple survey geometry where 
the scales considered are much smaller than the extent of the 
survey. 
The 



bias formal i sm (e.g. Huterer et al. l2006t 



Amara & Refr lgierl 125(31 iJoachimi & Schneider! l2009h allows 



us to compute the bias on the intrinsic alignment parameters due 
to the residual GG signal in the transformed power spectra via 



(33) 



ZZ n G^ C ° V ^<W' n GlW) 



e i.j 



dp Y 



Note that, contrary to works focusing on cosmic shear analy- 
ses and treating intrinsic alignments as the systematic, we use 
the GI contribution in d30b and insert the transformed GG signal 
into ( f33l > such that it plays the role of a systematic. We make the 
assumption that, given z, < zj, the galaxy redshift distribution 
p(''(z) entering (0 is sufficiently compact that we can take the 
term [(1 + Z/)/(l + z p i v )] r out of the comoving distance integra- 
tion. Then both parameter derivatives needed for < f30l > and ( 1331 ) 
are readily calculated analytically. 

Since at this point we merely seek to demonstrate the con- 
cept of boosting, we limit the set of U^'\{) entering (f30b to those 
bins i which fulfil z, £ [0.4; 1.4]. This ensures that the approxi- 
mation ( 1221 can be used throughout and that it is straightforward 
to assign o- opt . Besides, we avoid issues at high z; with non-zero 
G w (y(z max )), which could possibly violate the basic condition 
G w (^hor) = 0, see Sect. 12. 31 We determine cr opt by computing 
the diagnostic £ given by ( l25l l for all z, and devising simple, 
piecewise linear formulae which yield a cr opt in the regime of 
small f for every z,. For survey P2 (standard photo-z) we use 
cr opt = 0.02 Zi + 0.035. The two other surveys have the same red- 
shift binning and hence identical f. We set cr opt = 0.13 z/ -0.014 
forz, < 1 and <x opt = -0.057 z, + 0.173 forz/ > 1 in these cases. 

Now we are in the position to compute the boosting transfor- 
mation for power spectra with z,- 6 [0.4; 1.4]. By means of < f3Qb 
and (f33j we obtain statistical and systematic error estimates for 
both intrinsic alignment parameters for all three survey models, 
summarised in Table [3] When varying both parameters, we find 
marginalised 1 cr errors of approximately 2.9 for A and 7.4 for 
y in case of survey S. The two surveys with photometric red- 
shift data produce errors around 0.7 on A and of the order 2 for 
y. As expected, the bias due to the remaining cosmic shear sig- 
nal is negligible in the case of the spectroscopic survey S and 
clearly subdominant in the case of survey PI. Even for the stan- 
dard photo-z setup P2 biases remain within the statistical 1 cr 
errors, reaching up to |£> S ys/c s tatl ~ 1/4 for A. 

In Fig. [5] the corresponding 1 cr confidence contours in the 
parameter plane A - y are given for the three survey models. As 
we have chosen a pivot redshift which is below the minimum 
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Fig. 5. Constraints on the free parameters of the GI model. 
Shown are the 1 <x confidence contours for the different survey 
models used. The solid ellipse with its centre indicated by a filled 
circle corresponds to the spectroscopic survey S, the long-dashed 
ellipse and open circle to survey PI, and the short-dashed ellipse 
with diamond to survey P2. Note that the centres of the contours 
are offset due to the bias by the residual GG signal. The cross 
marks the fiducial values of A — 1 and y — 0. 

redshift of GI signals that enter the analysis, a positive y leads 
to an increase in the amplitude of the GI model, which can be 
compensated by a smaller A. Hence, A and y are anti-correlated, 
leading to the degeneracy as indicated by the error ellipses. The 
bias acts mainly on A because a residual GG signal will to zeroth 
order affect the overall amplitude of the signal. In all three cases 
the 1 <x contours comfortably enclose the fiducial, true parameter 
values. 

Due to the low number density of galaxies, survey S is 
clearly not competitive in constraints on intrinsic alignment 
properties. The results from the two other surveys are not ca- 
pable of pinning down the intrinsic alignment model with high 
precision, but their bounds are comparable to current constraints 
by analyses of spectroscopic meas urements of galaxy num - 
ber density-shape cross-correlations (Mandelbaum et al. 2010). 
Note that the weights used for this analysis may still have con- 
siderable room for optimisation, and that we only used a limited 
range of z, . 

Table [3] also lists the resulting errors when only A is varied 
and no additional redshift dependence of the intrinsic alignment 
model is assumed. Constraints improve significantly when lift- 
ing the degeneracy with y such that A is determined to better 
than ±0.3 (1 cr) for the survey models with photometric redshifts 
while constraints by survey S are about three times weaker. The 
bias is still negligible for the spectroscopic survey model, and 
clearly subdominant for survey PI (cr ph = 0.03). The residual 
systematic affects the error budget noticeably for the analysis of 
survey P2 (cr ph = 0.05) with |b sys /cr stat | » 63 %. Again, optimisa- 
tion of the boosting procedure may further decrease the residual 
cosmic shear signal well below the statistical 1 cr-limit. 

The errors for the good photo-z and in particular the spec- 
troscopic survey models are dominated by shape noise due to 



the low number density of galaxies in each tomographic galaxy 
sample, apart from only the smallest angular frequencies. As can 
be seen from (1321 1. the errors scale inversely with the total num- 
ber of galaxies in the survey if cosmic variance is negligible. 
Thus, if in the future larger number densities of galaxies with 
highly accurate photometric redshifts than assumed in this work 
are attainable, the constraints on GI correlations via the boosting 
technique will improve accordingly. If we re-run the analysis for 
survey S with the galaxy number density assumed for survey PI, 
i.e. a factor of 10 higher, all the statistical errors indeed decrease 
by almost an order of magnitude. 

6. Relation to GI nulling 

If one is able to extract the GI signal from cosmic shear data, the 
question arises whether this could also be used to remove the GI 
contamination from the data and thus make cosmic shear analy- 
ses robust against biases due to intrinsic alignments. Intuitively, 
one can simply subtract an isolated GI signal from the origi- 
nal measures, and indeed we are going to devise such a pro- 
cedure. Afterwards we will again propose a simple, paramet- 
ric weight function to construct a boosting method, whose out- 
come will then be used to eliminate the GI signal. These steps 
are not optimised and merely serve to demonstrate the link be- 
tween GI boosting and its removal, as well as to compare the 
performance of the la t ter to the st andard nulling technique of 
iJoachimi & Schneider! (12008, 2009) in a simple scenario. 

6.1. Signal transformation 

As an alternative to the procedure in Sect. l2.2l one can choose 
the lower integration boundary in (O as Xmin - Xi- As is evident 
from (O, in this case only the first term of the transformed GI 
signal remains. Hence, it is likely that^ m ; n = produces a larger 
amplitude of the modified GI power spectrum, but Xmm - Xi 
results in a cleaner signal insofar as it contains only contribu- 
tions from intrinsic alignments generated by matter at distance 
Xi- Consequently, we are going to use the latter choice of^ m ; n for 
constructing a method to remove the GI signal at^,. The trans- 
formed lensing signal for;^ m i n = Xi is derived in analogy to ( ITTb 
and reads 

o« = *I * i,0c 4-f) ,34) 

Now suppose we are able to construct a boosting technique 
with a sig nificant signal n® (0 while n®j(0 * 0. Noting again 
that the remaining first term in ((9]) is a rescaled version of the 
original GI signal (0, we define a further set of power spectra 

e*w - - fu n oLw with fij = . (35) 

and likewise for the individual GG and GI signals. This defi- 
nition holds for all i < j. The auto-correlations Q ( "\£) would 
simply correspond to the original auto-correlation power spec- 
tra P iu \(). As we are still working in the approximation of very 
narrow redshift bins, auto-correlations are hardly affected by GI 
correlations at all. In practice, auto-correlations are likely to be 
excluded or specially treated anyway due to the presence of in- 
trinsic ellipticity correlations, see the discussion in Sect.|7] 
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Assuming that the GI boosting works effectively, 



0, so that one expects that Qqq({) ~ P^Mf), i.e. the transformed 
cosmic shear signal is close to the original GG term. Switching 
to the notation of narrow redshift bins again, we find for the 
transformed GI signal 

Qci(Xi.Xj>0 = ~ 2 siXhXi) Psi\—,Xi\ (36) 

Xi \Xi I 

-fa - 2 G (xd Psi —,Xi = , 

2c 2 Xi XXi I 



where we have inserted ((TJ and the first term of ©, and made use 
of the transition g^(xd — * g(Xj>Xi)> see ©■ As a consequence, 

Qo&W * ^mW-^^W * *pgW- Hence ' if we can devise 
an effective boosting technique using Xmia = Xi< we immediately 

have a means of GI removal at our disposal via ( 1331 ). 

Note that the standard n ulling technique as presented in 

Joac himi & Schneider! ((2008) also makes use of the definition 

® with 

Xmin — Xi- The central condition in their approach is 
recovered in our formalism by requiring G w (y,) = 0, which 
eliminates the GI signal under the same assumption of narrow 
redshift bins, see For practical purposes we also switch to 
the discretised form of the signal transformation ( TT2b . using now 



6.2. Construction of weights 

We begin by developing again a boosting technique, now for the 
changed condition Xmin = Xi- Due to the associated change in the 
lower boundary of integration in ([8]), the condition to remove the 
GG signal is altered as well. Keeping the same approximations 
as used to derive ( |2TT i. we now obtain from ( 1341 ) 



X 



X 



"Xi /"".Mm 

H) J X , \ Xl\ Xi 
Xi 
2 



(37) 



Jo Jxi 

(Ml - M 2 y) = , 

where we executed the integration over^- and defined 

M,= A X B^(x) x 1 '* 1 ; A* =1,2. (38) 

Inserting ([T3l into the foregoing definition and integrating by 
parts, one arrives at the useful relations 



Mi = G^ixd-xi -s— to); 



(39) 



M 2 = — - — (xd 

dx 



When these are plugged into (l37l i. we obtain a condition which 
is the equivalent of (fJJJi, i.e. which ensures the suppression of 
the GG signal in the transformed power spectra (0, 



^-(Xi) = ^-GU(xd. 
ox 2Xi 



(40) 



In contrast to (fJJJl, which is an integral condition on G (,) (x) and 
in its discrete form involves all sampling points between imm 
and Zi, (l40b is local and even contains a derivative of G (l) (x). 
Hence, we suspect that ( |40T > is less robust against the inevitable 
discretisation of the weight function. As a cross-check for the 



accuracy of ( 1371 ). and equivalently d40l ). we define the additional 
diagnostic 



n = 



M!-M 2 j 



(41) 



^B (i) (x(Zj))x'(Zj)Azjil 



X(zi) 



where the integrals ( 1381 ) were transformed to redshift and discre- 
tised in analogy to $Y2\ . 

Moreover, ( f4Qb hinders us to impose the condition 
dG^ I dx\ Xi = again, which boosted the GI term, see (|9). We 
define 



G%(x) = N exp }[\ h) 



(42) 



which has one free parameter less than ( fTol ). To avoid any confu- 
sion with foregoing usage, we will add a sub- or superscript Q to 
indicate quantities which are used in this section for devising a 
nulling procedure. The condition (l40l readily implies b = ^,/3. 
As long as cr/xi <K 1, (l42l has an extremum in the vicinity of xu 
located at 



_ 2 1 

z^extr — ^' 




(43) 



Therefore Gq(x) as defined in d42l) should nonetheless boost 
the tranformed GI signal fairly well. In complete analogy to the 
derivation in Sect. [3] one obtains the weight function 



D (0/ ^ Kr 2 x f (x~XiY 

B qM = N— 2 exp — 



(44) 



x 20r 



X±, (X-Xi) 



a 7 
3^ + -x, 



The normalisation N is again given by ( fT9l ). As before, this 
weight function still has one free parameter cr which will be 
used to optimise the representations of the continuous functions 
d42l and d44b by the discrete set of sampling points entering ([SJ. 
The weights derived from d44l) yield GTboosted power spectra 
n ( o'(^) via ©i which in turn produce GTnulled measures via 

G3. 



6.3. Nulled signals 

Again we study a set of diagnostics as a function of cr, to identify 
regimes of cr z where the GI nulling performs well. In Fig. [6] we 
plot ( as defined in (l25T l. rj which assesses how ( |40T > is affected 
by the discretisation, and \Gq(Xi)\ as an indicator of the boosting 
of the GI signal in the IE; (^Q f° r the spectroscopic survey S at 

Zi = 0.53. Furthermore we show the GI over GG ratio rfL, which 
is given by (l28T l when replacing X by the nulled power spectra 
(l33T >. Note that small values of are indicative of an effective 
removal of the GI signal. 

One might expect that |Gg (xi)\ is largest for small cr- be- 
cause G ( q(x) is sharply peaked with a large maximum value. 
However, this effect is counteracted by the normalisation of the 



7 Note that since we have normalised Gq(x), \Gq ,(Xi)\ ls a meaningful 
measure of the size of Gq (x) at^,-, relative to its overall amplitude. 
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Fig. 6. Determination of cr op t for the spectroscopic survey S at 
Zi = 0.53. Top panel: Diagnostics ( (dotted line), 77 (solid line), 
and \Gq (xd\ (dashed line) as a function of cr z . Larger cr z yield 

the desired small values of rj and larger \Gq (xd\, m agreement 
with the minimum of f which is found at 0% ss 0.28. Note that 
we have rescaled ( for convenience. Bottom panel: GI over GG 



ratio r% as a function of <x z . The diagnostics indeed hint at a 



regime of cr z where r9. is smallest 



weight function. Large values of o~ z cause G^(x) to be smoother, 
i.e. to have smaller curvature. Due to (PT3l the amplitude of 
Bq(x) would thus decrease for fixed normalisation. Since we 

normalise B ( q(x) according to (T% for every cr z individually, 
large <x, yield a higher normalisation relative to small <x,, imply- 
ing also larger values of Gq (x)- Hence, one observes an increase 
in \Gq (Xi)\ as a function of o~ z . 

The diagnostic 77 has relatively large values for strongly 
peaked G ( q(x) and decreases slowly for larger cr z . A small 
change in the weight function Bq (x) due to the discretisation 
can induce significant changes in Gq(x), and its slope close 
to Xu which are the stronger the more sharply peaked Gq(x) 
is. Therefore d40b is more difficult to fulfil at small 0%. Both 
|Gg (xd\ and j] prefer larger cr z , in agreement with which we 
thus continue to use for the determination of cr opt . As £ clearly 
disfavours o~ z > 0.3, we choose as the optimum the minimum of 
( at about 0.28. Considering the lower panel of Fig. [6] this value 
is in very good agreement with small and hence close to optimal 
values of r9 v Generally, we find that cr opt is considerably larger 
for this approach, compared to the variant analysed in Sect.[5] 

With this finding at hand, we can compute Gl-boosted power 
spectra according to (O, and from these sets of nulled power 
spectra via (l35l l. results for both being shown in Fig. [7] The GI 
term is significantly less boosted than in the version studied in 
Sect. l5.1l with rci less than 10 (see Table[2]for comparison). Still, 
the intrinsic alignment suppression works excellently with rQ. < 
5xl0~ 4 for all background redshift bins and angular frequencies. 
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Fig. 7. Nulling performance for the spectroscopic survey S at 
Zi = 0.53. Top panel: Transformed GG (solid curve) and GI 
(dotted curve) power spectra, computed according to ([8J. Centre 
panel: GI over GG ratio as a function of angular frequency 
for the same parameters as above. The background redshift bin j 
is set to j = i+ 1 (solid black curve), j = (i+N z )/2 ('mid'; dotted 
black curve), and j = N z ('far'; solid grey curve). Lower panel: 
Relative deviation Agg as a function of angular frequency for the 
same parameters as above. The coding of the curves is identical 
to the foregoing case. 

In Table [4] values of for other z, and in addition for the 
good photo-z survey PI are listed. The down weighting of the 
GI signal quickly deteriorates with the increase in photometric 
redshift uncertainty, being more than two orders of magnitude 
larger for survey PL For the standard photo-z case we find that 
the boosting as implemented in this section is ineffective, so that 
we do not consider it here. As shown in Sect. 16. 11 under idealistic 
circumstances one expects the signal in the nulled power spec- 
tra Q ( 'i\{) to be close to the one in the original power spectra 
P^^(€). Hence, we calculate the quantity 



Agg = 



6gg(D 



1 . 



(45) 



which is also given in Fig.|7]and Table|4] The deviation from the 
original signal is at the per cent level for close foreground and 
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Table 4. Summary of the nulling performance for two survey models and different values of zi- Given are the median values of the 
GI over GG ratio of the power spectra Q^' J \(), r® v and the relative deviation of Q { '^{€) from the original power spectra, Agg, for 
the correlation of adjacent bins ('/ + 1'), the cross-correlation between bin i and the background bin with index (i + N z )/2 ('mid'), 
and the cross-correlation between bin i and the most distant bin with index N z at z < 2 ('far'). In addition cr opt as determined from £ 
is listed for every case considered. 



survey 


Zi 


^opt 


#&(»+ 1) 


r^(mid) 


rgj(far) 


AqgO'+I) 


A GG (mid) 


A GG (far) 


S 
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background redshift bins with z; < zj, and increases to about 
20 % if bins i and j are far apart, irrespective of the photometric 
redshift quality. 

6.4. Information content 

How does the nulling technique as outlined above perform in 
comparison with the standard nulling approach? For a very dense 
binning in redshift both methods evidently remove the GI con- 
tamination of the cosmic shear si gnal to h i gh ac curacy, see 
for instance the recen t findin gs bv IShi et all (12010). However, 
iJoachimi & Schneider! ((2009) have shown that, even in idealistic 
situations comparable to our spectroscopic survey, a substantial 
loss of cosmological information is inherent to standard nulling. 
We assess the information content in both nulling approaches in 
a simple case study. 

We restrict ourselves to the spectroscopic survey model S 
and consider again only zi € [0.4; 1.4], for the same reasons 
as discussed in Sect. 15. 21 Again, we compute £ for all z; to 
find a simple prescription for cr opt ; in this case we use cr opt = 
-0.131 z, + 0.346. The information content is quantified in terms 
of the cumulative signal-to-noise (S/N), defined as 

^ZZ cov- 1 (e«g(o. eM flj»o> > (4€> 

l j>i, l>k 

where the covariance of the nulled power spectra can be derived 
from ( 1331 ). 

Cov(Q m (0,Q ikl \()) = Co\(pfo\e), P <kl) (ij) (47) 
- Cov (n®0), P w (€)) - fa Cov(p m (.€), Ylf(€)) 

+ fnfa cov (n®(0, ng>(0) . 

The S/N for data sets of original pow er spectra P <lJ ^({) and 
of nul led power spectra obtained via the Joachimi & Sch neider! 
(2009) formalism are calculated in analogy to d46*b . For this setup 
it is safe to assume that d32l has only contributions from shape 
noise. Even with this simplification, the inversion of the covari- 
ance is computationally expensive for the total of 65 bins be- 
tween Zi = 0.4 and z, = 1 .4. Thus we include by default only 
tomographic measures for every fifth bin i, but all j > i, in the 
S/N. The absolute value of the S/N depends of course on how 
many power spectra are incorporated, but we are only interested 
in the ratio of S/N for the nulled datasets over the set of original 
power spectra. 

Note that for every z; one can make use of N z - i power 
spectra Q (, ^ > (€). The very same number of modes is available in 
the standard nulling approach although one mode is discarded to 



Table 5. Ratio of cumulative signal-to-noise of the nulled set of 
power spectra over the original set of power spectra (SNR). The 
results for the nulling m ethod devised in this work an d the stan- 
dard nulling technique (Joachimi & Schneider! 12009) are com- 
pared. The default SNR shown is computed for a step size in 
foreground redshift bin index i of 5 and in the redshift range 
[0.4; 1.4]. In addition the moduli of the fractional deviation of 
the SNR from the default when using a step size of 3, i.e. every 
third bin i (third column), a step size of 7 (fourth column), and 
the default step size but a different redshift range [0.5; 1.3] (fifth 
column) are given. 

nulling type SNR step 3 step 7 z e [0.5; 1.3] 
this work 0.179 0.76% 2.75% 1.96% 
standard 0.163 0.59% 2.32% 1.45% 



perfo rm the actual nulling (for details see Joachi mi & Schneider! 
120091 Transformed auto-correlation power spectra with i - j do 
not enter the S/N, but by construction the P^"\{) do contribute 
to all Q^{€) via the YTq(€), whereas in standard nulling auto- 
correlations are completely discarded. However, due to the dense 
redshift binning, we expect the amount of independent informa- 
tion contained in auto-correlation power spectra to be small. 

We have given the resulting ratios of the S/N for the nulled 
data set over the S/N for the original one in Table|5] The consid- 
erable loss of information can be confirmed, the S/N for both 
nulling methods yielding less than 20 % of the original S/N. 
We find that these numbers are very robust against changes in 
the number and values of redshift bins i included in the S/N by 
varying the size of steps in bin numbers i and the range of red- 
shifts considered. It is quite remarkable that the ratios for both 
nulling methods are very similar. The slightly bigger number for 
the nulling as devised in this work could be related to the inclu- 
sion of auto-correlation power spectra, but is not very significant 
anyway. 

In the standard nulling case the information loss is caused 
by discarding part of the signal, namely one mode per bin i 
whereas the variant suggested here features a signal that devi- 
ates by at most about 20 % from the untransformed one. In the 
latter case the loss is caused by an increase in the covariance due 
to the subtraction of signals in ( T35l l. We conjecture at this point 
that the agreement in the amount of information lost, in spite of 
the largely different mechanisms of the two methods, hints at a 
fundamental limit of how far GI and GG signals can be distin- 
guished by only relying on the redshift dependence of the two 
contributions. 
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7. Conclusions 

In this paper we presented a method which extracts shear- 
ellipticity correlations (the GI signal) from a tomographic 
cosmic-shear data set. The approach relies neither on models of 
intrinsic alignments nor on knowledge of the cosmological pa- 
rameters that characterise the cosmic shear (GG) signal, making 
only use of the typical and well-understood redshift dependen- 
cies of both the GI and GG term. We derived constraints which 
a linear transformation of second-order cosmic shear measures 
has to fulfil in order to boost the GI signal and simultaneously 
suppress the lensing contribution. We studied in depth a partic- 
ular parametrisation of the weights entering this transformation 
and analysed the performance of the resulting GI boosting tech- 
nique for three representative survey models. 

Applying the GI boosting to future all-sky cosmic shear sur- 
veys, it should be possible to isolate the GI signal with subdom- 
inant biases due to a residual GG term, and with constraints that 
are comparable to current res ults from indirect measur ements 
of shear-ellipticity correlations dMandel baum et al. 2010). If one 
restricts the analysis to galaxies with photometric redshift infor- 
mation of good quality, i.e. a redshift scatter of not more than 
o"ph(l + z) with <r ph = 0.03, one can achieve 1 cr-errors on the GI 
signal amplitude A in the parametrisation of (l29l > of better than 
0.2 when varying only the amplitude, and a marginalised error 
of approximately 0.7 when fitting an additional redshift depen- 
dence. 

Using all galaxies from a survey fulfilling <x ph < 0.05, the 
statistical constraints degrade only marginally but the parameter 
bias due to the residual GG contribution can attain more sig- 
nificant values of up to £> sys /<x stat < 2/3. We also considered 
a survey with high-quality photometric or spectroscopic red- 
shifts. However, the expected low number density of galaxies 
of n„ — 1 arcmin~ 2 , even for future surveys, does not permit us 
to place competitive constraints on intrinsic alignment models. 
In this case of highly accurate redshift information the residual 
bias on parameters is negligible. 

Although we have modelled scatter in photometric redshifts 
for our investigations, we did not consider other effects affect- 
ing the accuracy of redshift information, such as an error in the 
median of the galaxy redshift distributions or catastrophic fail- 
ures in the determination of photometric redshifts. As several 
studies of in trinsic alignment removal techniques have demon- 
strated (e.g.lBridle & Kind 120071; Uoachimi & Schneider! f2009t 
iJoach imi & Bridle 2010), the ability to separate the GI from the 
GG signal depends vitally on these parameters characterising the 
accuracy of and knowledge about redshifts. The same can be ex- 
pected for the GI boosting technique, possibly to an even larger 
extent since in this case one attempts to suppress the originally 
strongest contribution to ellipticity correlations, the GG signal. 
Hence, we hypothesise that the requirements of future ambitious 
weak lensing surveys, like a negligible fraction of catastrophic 
failures and an error in the m ean of each re dshift distribution of 
not more than 0.002(1 + z) (lLaureiis et alJ l2009'). are both nec- 
essary and sufficient for a success of GI boosting. We leave a 
detailed assessment of the requirements on the quality of red- 
shift information to future work. 

Moreover, we did not yet include intrinsic ellipticity corre- 
lations (II) into our considerations. Since the II signal is gener- 
ated by physically close pairs of galaxies, it has a redshift de- 
pendence that is clearly distinct from t he GI and GG terms, and 
can thus be removed relati vely easily (King & Schne iderll2002l 
120031 iHevmans & Heavensll2003t Takada & Whitdl2004l In to- 
mographic cosmic shear data it mainly affects auto-correlations 



and cross-correlations of adjacent photometric redshift bins with 
significant overlap of their corresponding distributions of true 
redshifts. One of the aforementioned II removal techniques 
could precede the GI boosting, causing an increased shape noise 
contribution in particular in the auto-correlations due to the re- 
duced number of available galaxy pairs. Alternatively, the down- 
weighting of the II signal could also be readily incorporated into 
the boosting technique by introducing the additional condition 
<9 2 G® /dx 2 \ X i = 0, implying B^(xi) = and therefore a down- 
weighting of auto-correlations as well as cross-correlations of 
adjacent redshift distributions, see ©. 

Our findings still have the potential for significant improve- 
ment because we have only considered one specific parametrisa- 
tion of the weight function that governs the boosting transforma- 
tion. While this choice is intuitive and allows analytical progress, 
a more versatile approach could be to assume the weight func- 
tion B (,) 0f) as piecewise linear, with nodes placed at the median 
redshift of every galaxy redshift sample. The constraints on GI 
boosting and GG suppression could then be directly imposed on 
the discretised version of the boosting transformation, thereby 
fixing a subset of the values of B^ (x) at its nodes. The remain- 
ing freedom in the weight function could for instance be used 
to maximise the signal-to-noise of the expected transformed GI 
signal. 

We also constructed a method of GI removal, directly based 
on a slightly modified version of the GI boosting technique. In 
principle, we showed that if one is able to isolate the GI sig- 
nal via boosting, one can simply subtract a rescaled version of 
the GI term from the original cosmic shear measures to elimi- 
nate the intrinsic alignment systematic. We find that the residual 
contamination of the cosmic shear signal by GI correlations is 
indeed small, and that the cumulative signal-to-noise of the thus 
treated cosmic shear signal decreases by about a factor of 6. This 
value is remarkably close t o the result for the standard GI nullin g 
technique as introduced by Joachi mi & Schneider! ( 120 08. 2009), 
in spite of the differing approaches. The underlying reason for 
this agreement may be due to a fundamental limit in the ability 
to separate GI and GG signals relying only on the dependence on 
redshift, which is worth to be addressed in future investigations. 
Of course, such a limit would also imply a maximum accuracy 
with which parameters of intrinsic alignments can be constrained 
via GI boosting. 

Like the method devised in this work, the standard nulling 
technique is also a purely geometrical method. Hence, a com- 
bined application of GI boosting and nulling to a cosmic shear 
data set would still be based on a minimum of assumptions about 
the actual forms of signals or the values of model parameters. 
For instance one could use an initial analysis based on nulling to 
yield robust estimates of the cosmic shear signal and the corre- 
sponding cosmological model. This could then be used to con- 
struct weights for the GI boosting transformation such that even 
in the case of standard photometric redshift quality (which we 
assumed to be cr ph = 0.05 in this paper) the bias due to the resid- 
ual GG signal would be negligible, thereby enabling an equally 
robust estimate of the GI signal. 

Ultimately, the cosmic shear analysis, the treatment of intrin- 
sic alignments, and the inclusion of addition al information from 
galaxy number density correlations (as in Mandelbau m et al.l 
2006: lHirataeTaIll2007l: lMandelb aum et alj|2010b will all be ef- 
ficiently c ombined into a s imult aneous analysis of the form pre- 
sented in Bernstein (120091) and IJoachimi & Brid le (2010), pro- 
vided one can summon the computational power. Yet the model- 
independent, direct, and robust boosting technique, as well as 
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nulling and the combination of the two, will prove useful e.g. 
to provide reliable priors on the large set of parameters entering 
the integrative approaches and in addition serve as a valuable 
consistency check in cosmic shear analyses. 
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